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SN2014J is the closest supernova of type la that occured in the last 40 years. This provides an opportunity for unprece¬ 
dented observational detail and coverage in many astronomical bands, which will help to better understand the still un¬ 
known astrophysics of these supernovae. For the hrst time, such an event occurs sufficiently nearby so that also gamma 
rays are able to contribute to such investigations. This is important, as the primary source of the supernova light is the 
radioactive energy from about 0.5 M© of ^®Ni produced in the explosion, and the gamma rays associated with this de¬ 
cay make the supernova shine for months. The INTEGRAL gamma-ray observatory of ESA has followed the supernova 
emission for almost 5 months. The characteristic gamma ray lines from the ^®Ni decay chain through ^®Co to ^®Ee have 
been measured. We discuss these observations, and the implications of the measured gamma-ray line characteristics as 
they evolve. 
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1 Introduction 

Supernovae of type la are commonly agreed to originate 
from the thermonuclear explosion of a CO white dwarf 
star in a binary system ( [Hoyle & Fowler||196Q| ). Accretion 
of material from the companion star onto the white dwarf 
eventually leads to nuclear ignition, and the nuclear energy 
release from carbon fusion ignited in the central degener¬ 
ate region of the white dwarf occurs so fast that the white 
dwarf cannot adjust its structure by expansion, and rather 
is disrupted (see, e.g., [Hillebrandt et~ar]|2013| for a recent 
review). Different paths have been proposed of how the su¬ 


pernova could be initiated from binary interaction (e.g. Pier-j 
santi et al.|2014| ): (1) The mass of the white dwarf may be 
increased from accretion of material from the companion 
star to and above the Chandrasekhar mass limit of stabil¬ 
ity for degenerate matter; (2) a major disturbing event may 
occur on the white dwarf surface and cause the white dwarf 
interior to become unstable towards runaway nuclear carbon 
fusion, e.g. accretion of a major amount of mass in form of 
a colliding body or a material cloud, or a nuclear explosion 
of accreted helium on the surface. 

Nuclear fusion at high densities processes the white 
dwarf material to iron group nuclei, which are the most sta¬ 
ble configuration of nuclear matter. ^^Ni is a likely prod¬ 
uct of such explosive supernova nucleosynthesis under typ¬ 
ical central densities of white dwarfs of 10^ g cm“^ or 
above and temperatures of a few GK ( jNomoto et al.|1997] ). 
The nuclear flame in principle propagates through heat con¬ 
duction, but is strongly accelerated through turbulent wrin¬ 


kling caused by instabilities, rushing through the star. The 
flame propagation thus is faster than any hydrodynamic ad¬ 
justment time scale, and an explosion is initiated, with lit¬ 
tle expansion of white dwarf material. At densities around 
10^ g cm“^, this defiagration may turn into an explosion, 
and nuclear burning then competes with expansion of the 
material, resulting in some outer parts of the white dwarf 
not being burnt towards iron group nuclei, but only to 
intermediate-mass elements, or even left unburnt as carbon 
and oxygen mainly ( jMazzali et al.pOOT] ). Typically, in a 
supernova of type la we expect that about 0.5 M© of nucle- 
osynthetic ^^Ni are embedded in about 0.5-0.9 M© of other 
material ( Mazzali et al.|2007}|Scalzo et al.|2014t|Stritzinger| 


jet al.|200^ 

^^Ni is unstable and decays first to ^^Co after r ^ 8 
days, then from ^^Co to ^^Fe at r ^ 111 days. Initially, the 
supernova is expected to be still opaque to even gamma rays 
at MeV energies, converting this radioactive energy into 
emission at lower energy photons ( [Clayton et al.j 1969][Hoe 


[flich et al.|1998 ). As the supernova expands, more and more 
of the ^^Ni gamma-rays from radioactive decay thus should 
be able to leave the source region where the decay occurs, 
and be observable with gamma-ray telescopes ( [Isern et aL 
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[2008[[Sim & Mazzali|2008[ ). In particular during the second 
decay stage producing gamma-rays at 846.77 and 1238.29 
keV, the supernova envelope is expected to become more 
and more transparent, so that radioactive decay and increas¬ 
ing transparency result in a maximum of gamma ray bright¬ 
ness at about 90-100 days after the explosion. As the super¬ 
nova unfolds, the rise and fall of the gamma-ray line inten¬ 
sity provides unique information on the morphology of the 
inner ejecta as imprinted by the explosion and the restruc- 
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R. Diehl: Gamma rays from SN2014J 


turing of the supernova as radioactive energy is deposited in 
its thus-shaped interior ( [Dessart et al.|2014[|Hoflich|2006| ). 
This has been the objective for gamma-ray astronomy since 
these prospects have been proposed by |Clayton et al.| ( p^69| ). 


SN2014J was discovered on January 22, 2014 (Fossey 


et al. [2014 1 in the nearby starburst galaxy M82 at 3.3 


Mpc distance ( [Foley et al.|2Q14| )as a type la explosion ( jCaoj 
et al.||2014] ). The explosion date appears to be 14 January, 


UT 14.75, with 0.2 to 0.3 d uncertainty ( jGoobar et al.|2015t 
Zheng et al. 2014| ). The supernova brightness maximum 


(blue band) was reached about 20 days after the explosion 
( jGoobar et al.||2014^ . Many studies from radio through 
infrared, optical and X-ray wavelengths were initiated, to 
study hints for its progenitor in pre-explosion data, hints for 
a companion star in emission details, and to follow the evo¬ 
lution of supernova light and spectra. 


2 INTEGRAL and SN2014J 


The INTEGRAL space gamma-ray observatory of ESA 
( [Winkler et al.||2003] ) observed SN2014J from end Jan¬ 
uary until end June, accumulating about 7 Msec of to¬ 
tal exposure; some other brief observations interrupted this 
SN2014J campaign ( [Kuulkers[[2Q14[ ), and hnahy visibility 
constraints ended SN2014J observations. 

The ESA International Gamma-Ray Astrophysics Lab¬ 
oratory ’INTEGRAL’ was launched into space in 2002, 


and carries two main telescopes, one called ’IBIS’ ( Uber- 
tini et al.|2003 ) and optimised for imaging, and one called 
’SEE ( Vedrenne et al. 2003) and optimised as a gamma- 
ray spectrometer. Both INTEGRAL main telescopes use the 
coded mask technique for imaging gamma-ray sources: a 
mask with occulting tungsten blocks and holes in the held of 
view of the gamma-ray camera casts a shadow from a celes¬ 
tial source onto the multi-element detector plane, which in 
the case of SPI consists of 19 Ge detectors, packed densely 
and with hexagonal shape. The SPI camera thus can resolve 
sources with a precision of 3 degrees, and its high-resolution 
Ge detectors obtain a spectrum at few keV resolution of ce¬ 
lestial gamma-ray sources. The coded-mask shadowgrams, 
however, have to be recognised above a large instrumental 
background caused by cosmic ray interactions in the space¬ 
craft and instrument materials. 


The data measured with SPI consist of energy-binned 
count spectra for the 15 of its 19 Ge detectors of the SPI 
telescope camera which were operational during our ob¬ 
servations of SN2014J. The campaign for SN2014J of the 
INTEGRAL mission involved orbit numbers 1380 to 1428, 
mostly dedicated to SN2014J, with several short interrup¬ 
tions for technical reasons or monitoring of other sources, 
and one major gap between 23 April and 27 May. These 
observations were planned to cover mainly the rising part 
of the expected gamma-ray line emission, combined with 
a longer exposure at late times, when gamma rays should 
not be absorbed any more by supernova material. Then the 


total amount of ^^Co would be measured as it decays. Ini¬ 
tially, we were in doubt as to where best invest the ’target of 
opportunity’ time, considering that with even 2 Ms of obser¬ 
vations on SN2011fe (distance about 6.5 Mpc) we had not 
achieved any hints for gamma-ray lines ( [Isern et al.|2Q13| ). 
But a few days after its discovery the supernova type and 
distance of about 3.5 Mpc were clear, and we decided to 
start observing on January 31, which was 16.3 days after 
the explosion ( [Kuulkers][2Q14| ). Recognising hints for the 
expected lines in quick-look data, we were able to convince 
the INTEGRAL user group and time allocation committee 
to continue monitoring of the supernova until visibility con¬ 
straints terminated this opportunity on 26 June, 164.0 after 
the explosion. 


INTEGRAL observations are typically made as 3000- 
second long pointings of the telescope towards a particular 
direction in the sky region of interest, moving the telescope 
axis by 2.1 degrees to shift the shadowgram of the source 
in the detector plane for the next 3000-second set, and so 
on, hnahy collecting exposure of the supernova in a regular 
pattern of telescope pointings of a 5 by 5 rectangle around 
the direction of SN2014J. The held of view of SPI is c^30 
degrees of opening angle, and always saw the source during 
those observations, at varying aspect angles. In total, after 
cleaning for data contaminated by solar hares or other irreg¬ 
ularities, we collected 1816 telescope pointings, each with 
15 detector spectra, hence 27240 spectra, each covering the 
20-2000 keV energy band at 0.5 keV bin width. 


Our analysis method is based on a comparison of mea¬ 
sured data to models in the complete set of count spectra as 
observed. For that, we must convert the expected shadow- 
grams for SN2014J in our set of pointings into the expected 
SN2014J counts using the SPI imaging response function, 
and we must develop a model for the large underlying in¬ 
strumental background. We then ht the intensity scaling fac¬ 
tor of the expected supernova contribution, plus a set of scal¬ 
ing factors per time for the instrumental background, to the 
set of measured spectra. 


Our model for the instrumental background is derived 
from a detailed spectroscopic assessment of the long term 
background and detector behaviour, in which we make use 
of the various characteristic instrumental lines as they reflect 
isotopes and their de-excitation gamma-rays, modulated by 
degradations of the spectral resolution of each of our detec¬ 
tors and their recoveries through periodic annealing opera¬ 
tions. Here we can account for the physical nature of instru¬ 
mental background lines and of detector-speciflc spectral re¬ 
sponses, and compared to earlier analyses we now combine 
data across a broader range of energy and time periods suit¬ 
ably to build a self-consistent description of instrumental 
background and its variations. Continuum and line contri¬ 
butions to the background are separately determined, detec¬ 
tor responses and their degradations are determined from a 
combination of spectral lines and their longterm behaviour, 
respecting each detectors speciflc response, but combining 
signals for statistical precision as possible. We And that 
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many aspects of the instrumental background can be un¬ 
derstood and constrained from modeling a broader spectral 
range of typically 100 keV width. This overcomes limita¬ 
tions of modeling each energy bin separately from all oth¬ 
ers, as done before. 


We have analysed the first part of data in particular to 
search for early appearance of characteristic lines from ^^Ni 
decay, as this would provide information of how close ^^Ni 
may have been to the surface of the supernova. ^^Ni de¬ 
cay occurs at r=8.8 days, and thus the first few weeks of 
data are suitable to address this question. The results of this 
analysis have been published in |Diehl et al. ( 2Q14| ). The full 
dataset then has been analysed to search for the appearance 
of the ^^Co decay lines as the supernova becomes more and 
more transparent. This decay at r=lll days occurs at the 
time scale at which the supernova reaches gamma-ray trans¬ 
parency, and the light curves and spectral shapes in the char¬ 
acteristic decay lines are the objective of this study, which 
have meanwhile also been published ( [Diehl et al.|2015| ). In 
the following two Sections, we summarise and discuss both 
these studies, as they were presented at the AG conference 
in Bamberg in September 2014. 


3 near the surface and its early 7 -rays 

SN2014J observations with INTEGRAL started about two 
decay lifetimes of ^^Ni after the supernova explosion. It 
was therefore a great surprise and unexpected to see the 
two strongest lines from ^^Ni decay, as shown in Fig[T] 
to 1^ Both lines appear nearly at their laboratory energy 
values, line shifts are constrained to below 2100 km s“^. 
Also line broadening is modest and constrained to below 
6000 km s“^, indicating that the ^^Ni near the surface does 
not expand with the higher velocities characteristic for outer 
supernova material, thus providing a hint for its possibly dif¬ 
ferent origin. 

This spectacular finding is puzzling in view of stan¬ 
dard defiagration-detonation models. We therefore scruti¬ 
nised our data and analysis methods, to investigate possi¬ 
ble systematic uncertainties, given that the statistical signif¬ 
icance of both lines is just at the 3cr level. We show in Fig.[^ 
that instrumental background lines present a major analy¬ 
sis challenge, but our modeling of instrumental background 
appears to properly account for them, within expectations. 
But instrumental background dominates the fiux uncertain¬ 
ties in each spectral bin of our result, and both statistical 
and systematic limitations may allow for other interpreta¬ 
tions. Testing this (see Fig. |^, we cannot exclude that we 
might have been fooled by statistical excursions and the 
truth could also be a broader line shape, or multiple line 
components; we just report the most likely interpretation, 
exploiting with confidence the spectral capabilities of our 
instrument, and properly accounting for Poissonian statisti¬ 
cal fluctuations as we fit our spectra. Finally, we point out 
that both lines at 158 and 812 keV are independent mea¬ 
surements, and are seen to fade in intensity consistent with 


Revolution 1380 



Fig.l Gamma-ray spectrum measured with 
SPI/INTEGRAL from SN2014J. The observed three- 
day interval around day 17.5 after the explosion shows the 
two main lines from ^^Ni decay. In deriving these spectra, 
we adopt the known position of SN2014J, and use the 
instrumental response and background model. Error bars 
are shown as Icr. The measured intensity corresponds to an 
initially-synthesised ^^Ni mass of 0.06 Mq. 



Fig. 2 Same as above, for the 812 keV line. We also il¬ 
lustrate our discrimination of sky and instrumental back¬ 
ground, showing the SN2014J spectrum against a scaled 
raw data spectrum dominated by instrumental background 
lines. Evidently, the line from SN2014J appears offset from 
the centroid of a strong background line, but may be af¬ 
fected by it, as shown by the two high data points found at 
the position of another strong background line. Note, how¬ 
ever, that the SN2014J line more-consistently follows the 
Gaussian response of SPI to a true gamma-ray line. 
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Fig. 3 Illustration of line position and shape significances. 
Here we re-analyse the spectrum shown above with differ¬ 
ent parameters for centroid and width, plotting the line sig¬ 
nificance for each such trial. Evidently, the narrow and un¬ 
shifted 158 keV line is the most probable result, but a lingle 
broad line as well as two other satellite lines may also be 
possible signals hidden in our data. 



15 20 25 30 35 40 45 

Day after Explosion 


Fig. 4 The fading of both lines from ^^Ni is consistent 
with the 8.8 day decay time, although statistical precision 
is poor, and observations should have better started earlier. 
We compare or measured intensity variation in the early pe- 
riod with severa l models from the set provided by |The & 


Burrows 


( 2014| ). Models with ^^Ni on or near the surface 
appear a better match than the canonical ’W7’ model, al¬ 
though this remains speculative. 


what would be expected for ^^Ni decay (r=8.8 days) (see 

Fig.g. 

How can we make sense of this surprising result? 
A single-degenerate Chandrasekhar-mass scenario appears 
unlikely to us: X-ray flux limits ( Nielsen et al.||2014 ) and 
pre-explosion images ( Kelly et al.||2014| ) exclude a super- 
soft progenitor. A sub-Chandrasekhar mass model with a 
He donor, or a merger of two white dwarfs, may seem bet¬ 
ter models for SN2014J, and could also explain our obser¬ 
vation of ^^Ni in the outer layer of the supernova more 
readily. Moreover, this He donor progenitor channel is fa¬ 
vored for this SN la from population-synthesis / supernova 
rate arguments ( [Ruiter et al.|2014| ). Also a classical double- 


to observer | 


I 

I 



-20,000 km s“^ 


Fig. 5 Sketch of a ejecta configuration compatible with 
our observations. Helium may have been accreted in a belt 
before the explosion, producing a ^^Ni-rich belt at the sur¬ 
face of the ejecta. The gamma-rays can escape early from 
the belt material, while the ^^Ni embedded more deeply and 
created by the main explosion (black) is still at high optical 
depths an invisible. A dashed arrow points to the observer, 
reflecting a face-on aspect within about 45° as required by 
a small if any shift of the observed line centroid. 


deton ation explosion scenario ( Fink et al.|2010 Moll et al. 
20141 is inconsistent with our observations: a ^^Ni shell 


engulfing the SN ejecta would be expected, resulting in 
broad, high-velocity gamma-ray emission lines, and more¬ 
over such an outer shell is expected to have an imprint on 
optical and infrared observables, which are probably not 
seen in SN2014J ( |Goobar et al.|2Q14bl|Telesco et al.|2Q15] ). 

A modified version of such model may, however, guide 
us towards what might have happened (see Fig. [^: If He 
would be accreted rapidly and form an equatorial accre¬ 
tion belt before it detonates, instead of accumulating in a 
shell, the kinematic constraints could be met, provided we 
observe the binary system essentially perpendicular to its 
orbital plane. Such an idea had been discussed frequently in 
the context of classic al novae ([Kippenhahn & Thomas|1978[ 


Law & Ritter||1983||Piro & Bildsten||2004| ). This might be 

consistent with the observed gamma-ray signal, and com¬ 
patible with optical observables. Our radiation transfer sim¬ 
ulations in UV/optic al/NIR (see [Diehl et al.|2014[ based on 
radiation transport, ( [Kerzendorf & Sim||2014| )) shows that 
the Ni-belt would not produce easily distinguishable fea¬ 
tures but result in normal SN la appearance, not only for a 
pole-on observer but also for an equatorial observer. In view 
of this, our interpretation of an externally-triggered explo¬ 
sion may be plausible, though speculative. Further observ¬ 
ables need to be checked against such a type of model. 
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4 How 7 -rays from ^®Co decay are revealed 


Clearly, in all variants of supernova type la models it is 
quite likely that major amounts of ^^Ni are produced, and 
also significant amounts of unbumt material are ejected, 
which occult gamma-rays from radioactive decay in the first 
months, as the supernova is spreading out. Unclear remains 
if ^^Ni is mixed throughout the supernova in the explosion, 
and hence some of it may already appear near the observ¬ 
able surface at early stages; this has been suggested by 3- 


dimen sional hydrodynamical simulations ( [Seitenzahl et al. 
2013| . Also unclear remains if ^^Ni production in the cen 


tral parts is reduced due to electron capture dominance at 
high densities, as suggested by fiattened line profiles of Co 
lines that can be seen in the mid infrared ( [Gerardy et al. 


2007 1 ). All occultation of gamma rays should decrease with 


time and be rather unimportant after about three months. 


Other INTEGRAL data analysis results had reported 
^^Co decay with its strongest lines at 847 and 1238 keV to 
appear as broad lines, quite in concordance with standard 
Chandrasekhar mass models ( [Churazov et al.|[2014| ). Our 
own analysis of time integrated data confirms this general 
picture (see Fig.[^. 

But as we exploit the full spectral resolution of our in¬ 
strument, we find that these broad lines seen at late epochs 
are not quite the same (though just occulted and at lower 
brightness) towards earlier epochs, see Fig. [7] They may, in 
fact, be composed of narrower emission lines which vary in 
intensity as the supernova unfolds. This could refiect emis¬ 
sion from a few ^^Co-rich plumes, embedded in the su¬ 
pernova and thus progressively revealed through structured 
non-radioactive ejecta with a more complex morphology. 


It is difficult to cut our observations into as many spec¬ 
tral and timing bins as we would like, to disentangle the 
signature from how ^^Co may be embeded in the super¬ 
nova. Attempting to maximise time resolution without im¬ 
posing a bias from a particular supernova model or simu¬ 
lation, we generated a light curve from our data alone, in 
the more-significant signal of the 847 keV line (see Fig.[^. 
In this Figure, we show for comparison also the four epochs 
where Figure|^shows the detailed spectra and their intensity 
variations seen more significantly. Clearly, this borders on 
what can be extracted given the uncertainty of our measure¬ 
ments; but our aim is to learn from the specific messages 
of gamma-ray data. Alternatively, if 3-dimensional hydro- 
dynamic models could be combined with current sophisti¬ 
cated radiation transport models, we could compare such 
more fiexible 3D model-predicted gamma-ray light curves 
to our observations. This is, however, beyond the scope of 
current simulations and analysis. Our analysis shows that 
significant variability characterises how gamma rays from 
^^Co decay appear, as time goes on and the supernova ejecta 
become gamma-ray transparent, clearly beyond the smooth 
light curves currently available from supernova model pre¬ 
dictions. More details can be found in [Diehl et al.| ( |2Q15] ). 



Fig. 6 SN2014J spectrum near the 847 keV line {above) 

and near the 1238 keV line (below) as expected from ^^Co 
decay. These spectra are determined in energy bins of width 
10 keV over the entire observing period; the source intensity 
is fitted at four independent epochs. For illustration, fitted 
Gaussians indicate the detection of broadened lines near the 
^^Co gamma-ray line energies. 

5 Summary and conclusions 

Supernova SN2014J is the first supernova of type la which 
is close enough for significant measurements of character¬ 
istic gamma rays from the ^^Ni decay chain. INTEGRAL 
has followed the gamma-ray emission for five months, and 
thus obtained the first clear detection of such characteristic 
gamma ray emission. This provides the first direct confir¬ 
mation that radioactive decay from ^^Ni through ^^Co is 
the energy source of supernova light, and that about 0.5 Mq 
of ^^Ni have been seen in gamma-rays from SN2014J. 
Gamma-ray data directly from ^^Ni radioactive decay have 
finally confirmed a key aspect of our understanding of type 
la supernovae. 

Upon a closer look, there are challenges and surprises, 
presented by those same gamma ray observations: Early 
time ^^Ni gamma-ray emission has been found, and is sur¬ 
prising in appearance. Its brightness suggests a major frac- 
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Fig. 7 SN2014J signal intensity variations for the 847 

keV line {center) and the 1238 keV line {right) as seen in 
the four epochs of our observations, and analyzed with 10 
keV energy bins. The 1238 keV fluxes have been scaled by 
the ^^Co decay branching ratio of 0.68 for equal-intensity 
appearance. Clear and signiflcant emission is seen in the 
lower energy band {left and center) through a dominating 
broad line attributed to 847 keV emission, the emission in 
the high-energy band in the 1238 keV line appears consis¬ 
tent and weaker, as expected from the branching ratio of 
0.68 {right). Fitted line details are discussed in the text. For 
the 847 keV line, in addition a high-spectral resolution anal¬ 
ysis is shown at 2 keV energy bin width {left), confirming 
the irregular, non-broad-Gaussian features in more detail. 


tion of near 10% of the total ^^Ni be visible 17 days after 
the explosion. It is unclear if this is ^^Ni contributed by a 
surface event, or if this happened to be a ^^Ni rich plume 
rising early from the inner supernova region to the surface. 
The later gamma-ray spectra and their ^^Co decay lines then 
show structure that may raise signiflcant doubts about a ho¬ 
mogeneous and smooth distribution of ^^Ni throughout the 
supernova. It remains a challenge how these ^^Co decay 
data can be reconciled with data from the receding photo¬ 
sphere at lower wavelength regimes, which altogether build 
a tomographic view of the angle-averaged morphology of 
the supernova. SN2014J provides a challenge to our under¬ 
standing of type la supernovae, on second glance. We will 
learn if that remains a specific and puzzling event, or if un- 
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Fig. 8 SN2014J signal intensity variations for the 847 

keV line in two different time resolutions. The 4-epoch re¬ 
sults are consistent with 11-epoch analysis, both showing 
an initial rise and late decline of ^^Co decay line intensity, 
with a maximum at 60-100 days after explosion. Shown are 
also se veral candidate source models from lThe & Burro^^ 
(2014 ), fitted in intensity and thus determining ^^Ni masses 
for each such model. The best-fitting model is shown as a 
continuous thick line. The ^^Ni mass has been derived from 
such model fits as 0.49±0.09 Mq (see piehl et al.|2015l for 
details) 


derlying processes linked to ^^Ni radioactive decay may not 
be apparent in data of stable ejecta. 
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